#include <unistd.h>
#include <stdio.h>
#include <math.h>
#include <mpi.h>
#include <stdlib.h>

#define DEBUG

int modulo(int, int);
void MPI_Op_Call(double *, double *, double *, int, MPI_Op);
int gthpc_Allreduce(void *, void *, int, MPI_Op op, MPI_Comm comm);
